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GIANT  COMETS,  EVOLUTION  AND  CIVILIZATION 


S.V.M.  Clube,  F.  Hoyle,  W.M.  Napier,  N.C.  Wickramasinghe 


Climatic  cycles,  ice  epochs,  mass  extinctions  and  other 
global  disturbances  may  be  attributed  to  episodes  of 
bombardment  by  giant  comets.  Such  bodies  now  in 
chaotic  orbits  beyond  Jupiter  present  the  most  serious 
current  celestial  hazard. 


Evidence  for  the  periodic  incidence  of  large-scale  terrestrial  phenomena 
[1]  has  grown  substantially  during  the  past  70  years.  The  periodic  value  ~26 
million  years,  moreover,  is  strikingly  close  to  the  interval  between  successive 
Solar  System  crossings  of  the  Galactic  plane  [2],  a  fact  which  suggests  that 
there  might  exist  a  predictable  celestial  mechanism  which  exerts  a  control 
over  such  fundamental  processes  as  the  evolution  of  life.  In  spite  of  its  sci¬ 
entific  and  human  importance  [3],  however,  it  is  only  during  the  last  quarter 
of  the  present  century  that  the  physical  aspects  of  the  galacto-terrestrial 
relationship  have  ceased  to  be  enigmatic. 

Thus,  whereas  until  the  mid-1970’s,  the  Oort  cloud  was  generally  thought 
of  as  a  stable,  primordial  system  essentially  not  influenced  by  its  Galactic 
environment,  it  now  appears  that  sequences  of  giant  comets  perturbed  into 
Earth-crossing  orbits  by  periodically  varying  tidal  forces  of  Galactic  origin 
are  the  probable  key  to  terrestrial  evolutionary  processes  [4].  It  turns  out, 
for  example,  that  the  ~26  Myr  period  of  the  late  Phanerozoic  pulsation 
(i.e.  since  about  260  million  years  BP)  is  indicative  of  the  Galaxy’s  longer 
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term  influence  on  the  Oort  cloud  while  the  current  (Pleistocene/Holocene) 
ice  epoch  and  its  ~0.1  Myr  alternation  between  glaciation  and  greenhouse 
are  broadly  indicative  of  a  current  more  localized  Galactic  influence  close  to 
the  plane  [5,  13].  The  current  arrival  rate  of  giant  comets  implies  that  they 
are  a  prime  celestial  hazard  to  civilization;  thus  it  appears  that  the  last  (i.e. 
latest  Pleistocene)  glacial  and  its  aftermath,  the  current  (Holocene)  inter¬ 
glacial,  may  be  associated  with  the  disintegration  and  decline  of  the  single, 
most  recent,  giant  comet  in  near-Earth  space  [6].  This  has  led  to  a  very 
natural  focus  now  on  the  lesser  but  still  significant  hazards  due  to  the  still 
disintegrating  debris  of  this  giant  comet,  producing  dark-  and  mini-ice  ages 
on  millennial  and  centennial  timescales  respectively  [3].  The  near-Earth  ar¬ 
rival  frequency  of  giant  comets,  ~  10~5  per  annum,  is  much  the  same  as  that 
of  ~l-kilometer  asteroid  impacts,  but  the  adverse  environmental  effects  of 
giant  comets,  although  comparable  in  severity,  axe  much  more  prolonged.  In 
the  present  study  we  consider  the  giant  comet  issue  from  both  astrophysical 
and  climatic  perspectives.  A  point  at  issue  is  the  greater  complexity  of  the 
cometary  as  opposed  to  the  asteroidal  hazard,  as  well  as  its  more  pervasive 
effects. 


Periodic  Extinctions  and  External  Forcing 

The  coincidences  in  age  between  tektites  and  Tertiary  extinction  events 
were  an  early  indication  that  extinctions  and  cometary  impacts  could  well 
be  related  [7].  An  extinction  mechanism  then  proposed  was  the  attenuation 
of  sunlight  through  stratospheric  dust  injection  from  the  meteor  stream  of  a 
large  comet,  not  least  for  the  Cretaceous-Tertiary  event  [8,  40].  Furthermore, 
it  did  seem  that  both  impact  cratering  and  global  tectonic  events  could  be 
understood  in  terms  of  bombardment  episodes  every  ~26  million  years  [9]. 
The  specific  Galactic  connection  only  emerged  however  with  the  realization 
that  global  terrestrial  disturbances  associated  with  bombardment  episodes 
could  be  both  regular  (periodic)  and  random  on  appropriate  timescales  pro¬ 
vided  large  comets  were  correspondingly  dislodged  from  the  Oort  cloud  as 
the  Solar  System  moved  through  the  Galaxy.  Conceived  as  a  general  theory 
of  catastrophic  evolution,  this  meant  that  major  extinctions  such  as  the  KT 
event  were  probably  due  to  ~  10-kilometer  (cometary)  asteroids  [10,  5].  Thus 
the  implied  galacto-terrestrial  relationship  was  based  on  prior  astrophysical 
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considerations  and  was  by  no  means  purely  a  response  to  the  discovery  of 
iridium  anomalies  [11]  or  other  extraterrestrial  markers  at  the  KT  boundary. 

Appropriate  geological  databases  have  now  improved  to  the  point  where 
the  periodicity  claims  can  be  subjected  to  quantitative  analysis  and  accu¬ 
rately  assessed.  Thus  in  a  1992  compilation  [12]  of  databases  of  the  past 
~260  Myr,  seventy  six  well-dated  major  geological  events  were  catalogued 
comprising  11  mass  extinction  peaks,  13  ocean  anoxic  and  black  shale  events, 
12  continental  flood  basalt  vulcanisms,  9  large  discontinuities  in  sea-floor 
spreading  rates,  17  dates  of  widespread  stratigraphic  sequence  boundaries 
(usually  indicating  major  periods  of  low  sea  level)  and  14  mountain  building 
events.  Power  spectrum  analysis  of  these  data  yields  evidence  of  an  overall 
~26  Myr  periodicity  coupled  with  a  recent  out-of-phase  surge  [12,  13].  The 
latter  commenced  some  5  Myr  BP,  merging  with  the  onset  of  the  Pleistocene 
glaciations  about  2.5  million  years  ago. 

In  Table  I  the  ‘best  fit’  periodicity  for  all  the  data  combined  is  compared 
with  the  individual  mass  extinction  peaks  from  the  compilation,  and  with  the 
periodicity  for  all  the  other  geological  phenomena  (cf.  Fig.  1).  In  construct¬ 
ing  the  peaks  the  recent  significant  surge  of  activity  (in  the  Pleistocene)  has 
been  excluded.  Clearly  the  mass  extinctions,  and  the  other  major  geological 
events,  fit  closely  to  the  same  ~26  Myr  periodicity,  separately  and  together. 
The  amplitude  of  the  periodicity  is  large,  the  whole  Earth  switching  on  and 
off  in  a  regular  cycle.  The  period  and  phase  are  remarkably  stable  against 
progressive  truncation  of  the  dataset  (Fig.  2). 

The  existence  of  a  periodicity  in  the  major  episodes  of  geological  activity 
strongly  constrains  the  possible  forcing  mechanisms.  The  end-Cretaceous 
extinctions  (~63  Myr),  the  Deccan  Trap  outpourings  (~65  Myr  BP)  and 
the  Chicxulub  impact  crater  (65.0±0.2  Myr)  are  closely  coincidental  in  time, 
but  also  fit  on  to  the  64  Myr  BP  peak  of  the  ‘Holmes  cycle’.  Unless  these 
coincidences  are  due  to  chance,  it  appears  that  the  periodicity  cannot  be  ex¬ 
plained  by  endogenous  forcing,  say  arising  from  mantle  convection  processes. 
Likewise  the  impact  of  stray  bodies  from  the  asteroid  or  Kuiper  belts  seems 
not  to  be  primarily  involved  -  notwithstanding  the  Spaceguard  project  [14] 
-  since  these  belts  are  not,  so  far  as  is  known,  subject  to  disturbances  with 
this  periodicity.  A  further  constraint  arises  from  the  fact  that  the  dynamical 
lifetime  of  an  object  in  a  near— Earth  orbit  is  about  30  million  years,  ade¬ 
quate  to  smear  out  any  periodicity  in  their  arrival  rate  unless  the  perturbing 
Earth-crossers  are  physically  short-lived.  Rocky  or  iron  asteroids  are  there¬ 
fore  excluded  as  prime  movers,  consistently  also  with  the  absence  of  bulk 
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meteoritic  material  at  the  KT  and  other  extinction  boundaries. 

Extraterrestrial  3He  of  probable  cometary  origin  has  recently  been  de¬ 
tected  in  a  North  Pacific  pelagic  clay  core  dating  back  over  70  Myr  [15].  The 
implied  flux  of  cometary  dust  is  highly  variable,  showing  strong  peaks  at  ~66 
Myr,  ~37  Myr  and  <  2  Myr  BP  (c/ Table  I).  There  is  thus  direct  evidence 
for  an  association  between  comet  dusting  and  geologically  noisy  epochs.  At 
higher  resolution,  the  dust  deposition  correlates  strongly  with  a  0.1  Myr  cy¬ 
cle  in  the  climatic  record  of  the  past  million  years  [15].  The  latter  seems  to 
imply  a  likely  enhanced  capture  probability  of  comets  as  the  ecliptic  crosses 
mid-latitude  Galactic  bands  [19]  and  where  comets  have  perihelion  and  nodal 
precession  rates  ( g ,  /)  which  match  the  ‘grand  cycles’  close  to  simple  multi¬ 
ples  kj  (j= 5, 6, 7, 8)  of  the  major  planetary  secular  resonances  ( gj,Sj )  [16],  eg 
when  \g\,  |/|  ~  3.25  arcsec/yr. 


Galactic  Encounters  and  the  Pleistocene  Surge 

Such  constraints  and  the  longer  term  periodicity  are  however  in  quanti¬ 
tative  agreement  with  the  perturbation  function  applied  to  the  long-period 
comet  system,  or  Oort  cloud.  This  comprises  some  1011  objects  which  orbit 
the  Sun  at  distances  of  up  to  60  000  astronomical  units  and  remain  only  just 
gravitationally  bound  to  it.  It  is  through  the  action  of  an  externally  induced 
tide  (or  differential  gravitational  field  with  respect  to  the  Sun)  that  a  flux  of 
comets  into  the  planetary  system  is  maintained  [17].  That  part  of  the  flux 
which  is  deflected  rapidly  and  terminally  into  sub-Jovian  space,  i.e.  through 
Jupiter-family  and  Halley-type  orbits  [4],  becomes  the  most  likely  to  interact 
with  the  Earth.  It  follows  that  the  comet-Earth  interaction  varies  overall  in 
direct  proportion  to  the  ambient  differential  field  affecting  highly  elongated, 
Oort  cloud  orbits.  Such  tides  arise  in  general  through  Oort  cloud-cum-Solar 
System  encounters  with  the  other  relatively  moving  elements  of  the  Galactic 
environment;  they  comprise  both  steady  and  impulsive  components  and  so 
produce  both  periodic  and  random  surges  in  the  giant  comet  flux  (Fig.  3). 
Depending  on  the  assumed  scale  height  of  dark  matter  in  the  Galactic  disc, 
the  periodic  modulations  may  easily  attain  amplitudes  of  3:1  or  more.  Con¬ 
sistency  with  the  observed  Holmes  cycle  is  obtained  if  the  half-period  of 
the  solar  vertical  oscillation  is  ~26  million  years.  This  corresponds  to  an 
in-plane  density  pz  ~  0.18  Mq  pc~3  for  the  ambient  Galactic  disc  [18],  im- 
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plying  the  presence  of  dark  matter  in  the  disc  with  an  extreme  Population  I 
distribution. 

The  long-period  comet  aphelia  currently  possess  an  underlying  isotropic 
distribution  deficient  in  objects  close  to  the  Galactic  poles  and  equator, 
broadly  in  accordance  with  the  Galactic  tide  [19]  complemented  by  individ¬ 
ual  (stellar  and  dark  matter)  perturbers.  In  addition,  however,  there  is  an 
excess  of  cometary  aphelia  around  the  Solar  antapex  direction  spread  along  a 
90°  segment  of  the  great  circle  which  includes  the  Galactic  centre  [20].  This 
accords  very  well  with  the  current  Solar  vector  and  the  observed  location 
of  the  nearest  large  mass  in  the  Galactic  environment,  namely  Gould’s  Belt 
(Fig.  4).  The  latter  appears  to  be  a  disintegrating  molecular  cloud  of  one  or 
two  million  solar  masses,  which  the  Solar  System  began  to  encounter  10  ±  2 
million  years  ago.  Allowing  for  the  infall  time  (~3-4  million  years)  of  dis¬ 
turbed  long-period  comets,  passage  of  the  Solar  System  through  the  Belt 
may  be  associated  with  the  recent  Pliocene  surge,  which  began  ~5  Myr  BP 
and  includes  three  orogenic  events,  two  sea  floor  spreading  discontinuities, 
five  evaporite  depositions  and  a  mass  extinction  ~1.5  Myr  BP  merging  with 
the  onset  of  the  current  ice  epoch  [13]. 


Comet  Splitting  and  Stratospheric  Dust 

The  observed  inverse  semi-major  axis  (i.e.  orbital  energy)  distributions 
of  long-  and  short-period  comets  have  long  been  understood  in  terms  of 
strong  planetary  encounters  in  association  with  hierarchical  splitting  (cf. 
Comet  P/Shoemaker-Levy  9),  the  fragments  being  increasingly  faded  and 
resilient  [21],  as  would  be  expected,  for  example,  were  the  largest  represen¬ 
tatives  (diameters  ^  100  kilometers,  say)  to  have  carbonaceous  chondritic 
cores  just  as  the  larger  parent  bodies  of  meteorites  (diameters  >  500  km, 
say)  are  expected  to  have  more  purely  chondritic  cores.  Thus,  to  the  ex¬ 
tent  that  the  long-period  comets  are  also  the  source  of  short-period  ones, 
classified  as  Jupiter-family  or  Halley-type  [22],  it  is  natural  that  these  deriva¬ 
tives  should  continue  to  split  and  are  represented  as  a  whole  by  populations 
in  excess  of  the  expected  unsplit  dynamical  transfers  which  include  large 
devolatilized  members  such  as  Chiron.  Splitting  is  a  commonly  observed 
process:  21  comets  have  been  recorded  as  doing  so  between  1846  and  1976 
and  three  since  1989  [23].  It  occurs  at  a  rate  of  one  or  more  events  per 
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comet  per  century  and  is  a  likely  major  mode  of  disintegration  and  evolu¬ 
tion.  The  precise  mechanism  is  not  known  except  for  those  comets  which 
pass  within  the  Roche  limit  of  a  planet  or  the  Sun;  and  the  emergence  of 
differing  constitutional  types  among  remnants  on  the  occasions  of  these  rare 
events  seems  to  indicate  compositions  which  are  heterogeneous  on  all  scales, 
including  some  which  cannot  withstand  typical  tidal  or  explosive  forces.  It  is 
likely  that  short-period  comets  make  close  passages  to  Jupiter  every  century 
or  so  and  if  one  in  a  thousand  of  these  is  a  giant  comet  (diameter  >  100 
kilometers),  splitting  into  thousands  of  fragments,  then  major  upsurges  in 
the  Jupiter-family  population  will  take  place  within  its  dynamical  lifetime, 
thereby  accounting  for  an  observed  excess  in  this  population  [24].  A  more  ex¬ 
tended  disintegration  pattern  can  be  expected  for  the  giant  comets  deflected 
into  sub-Jovian  space  as  they  will  tend  to  experience  close  encounters  with 
terrestrial  planets.  For  example  sub-Jovian  space  currently  contains  a  single 
broad  complex  of  (cometary)  asteroidal  and  meteoroidal  debris  known  as  the 
Taurids,  with  a  dynamical  lifetime  of  ~  10-2  Myr  [6]. 

In  general,  visible  comets  grow  tails,  losing  dust  and  volatiles,  when  their 
orbits  acquire  perihelia  of  q  <  2.5  AU,  and  many  are  associated  with  meteor 
streams  observed  when  the  Earth  intersects  their  orbital  tracks.  The  active 
comets  are  clearly  evanescent  bodies,  with  lifetimes  one  to  ten  millennia,  and 
so  represent  the  tip  of  an  iceberg,  with  a  large  undiscovered  population  of 
inert  bodies  belonging  to  the  Jupiter  and  Halley  families.  The  mass  distri¬ 
bution  of  comets  is  a  power  law  with  index  ~-1.7,  possibly  steepening  to 
~-2  at  the  high  mass  end.  The  cumulative  flux  of  large  long-period  comets 
is  given  by  F  ~  1  x  (d/5)-2  comets  AU-1t/r-1,  d  the  diameter  of  the  comet 
in  kilometers  [4].  There  is  no  securely  known  upper  limit  to  d,  but  sev¬ 
eral  historical  comets,  such  as  the  Great  Comet  of  1577  and  Comet  Sarabat 
of  1729,  appear  to  have  had  diameters  in  the  range  100-300  km  and  masses 
~  1021  — 1023  g.  Chiron,  in  a  chaotic  orbit  which  currently  lies  beyond  Saturn, 
may  be  180  ±  30  km  in  diameter.  A  giant  long-period  comet  (d  >  100  km) 
is  therefore  expected  to  cross  the  Earth’s  orbit  about  once  every  400  yr,  a 
Sarabat-sized  body  ( ^  300  km)  once  within  the  timescale  of  civilization.  The 
injection  rate  of  giant  comets  from  a  chaotic,  trans-Saturnian  orbit  into  a 
stable  Earth-crossing  one  is  of  order  10  Myr-1. 

A  comet  with  d  =  200  km  and  mass  M  =  1023  g  thrown  into  a  Taurid-like 
orbit  P(~3.3  yr,  eccentricity  ~0.85)  will  lose  M  ~  1018  gyr~l  due  to  out- 
gassing,  more  than  half  of  it  as  meteoroidal  dust  with  diameters  in  the  range 
0.01  ^m  to  several  mm  [25].  A  particle  size  distribution  n(a)da  oc  a-3-6  da, 


6 


a  >  0.01  mm,  is  indicated  by  the  Halley  data,  implying  a  significant  pile-up 
of  mass  in  the  submicron  size  range.  Outgassing  and  dust  production  will  not 
be  uniform  with  time:  the  overall  active  lifetime  of  the  comet  ~3000  yr  may 
be  interspersed  with  dormant  periods  when  the  surfaces  become  temporarily 
crusted.  During  its  active  lifetime  such  a  comet  could  generate  a  zodiacal 
cloud  of  mass  ~300  times  that  of  the  present  one. 

Solar  radiation  pressure  forces  have  an  important  role  in  separating  the 
various  grain  sizes  and  compositions  within  the  cloud.  The  critical  ratio  of 
radiation  pressure  force  to  gravity  below  which  particles  will  be  retained  in 
bound  orbits  depends  on  the  eccentricity  of  the  parent  cometary  orbit  as  well 
as  the  precise  point  on  the  orbit  at  which  the  particles  are  released.  For  a 
comet  in  a  circular  (or  near  circular)  orbit  with  radius  1  AU,  the  condition 
for  retention  is  P/G  <  0.5;  for  a  comet  in  an  elliptical  orbit  with  e  =  0.7 
(a  Jupiter  comet)  the  appropriate  condition  is  P/G  <  0.85  for  aphelion 
emission,  and  P/G  <  0.15  for  perihelion  emission  [39].  These  conditions 
are  marked  by  the  dotted  and  dashed  lines  in  Fig.  5.  The  several  curves 
show  the  ratio  P/G  from  Mie-type  computations  for  spherical  particles  of 
graphite,  iron  and  dielectric  compositions  (minerals  or  organics).  We  consider 
dielectric  grains  with  varying  degrees  of  porosity,  D10o,  D60,  D40  referring 
to  non-porous,  60%  volume-filled  and  40%  volume-filled  respectively.  The 
dielectric  material  in  these  calculations  is  taken  to  have  a  bulk  refractive 
index  m  —  1.5  and  density  2  g  cm~3.  Whilst  sub-micron  iron  and  graphite 
(metallic)  grains  are  seen  to  be  easily  lost  according  to  our  computed  P/G 
criteria,  organic  or  mineral  grains  of  all  sizes  which  are  <60%  volume-filled 
would  be  wholly  or  mostly  retained  in  bound  orbits,  depending  upon  the 
point  in  the  orbit  at  which  they  are  released. 

For  dust  particles  in  bound  orbits  the  main  loss  mechanism  from  the  zo¬ 
diacal  cloud  would  be  due  to  Poynting-Robertson  drag  forces.  The  Poynting- 
Robertson  timescale  at  1  AU  is  plotted  in  Fig.  6  for  each  of  the  grain  types 
considered.  We  note  that  iron  grains  of  radii  0.1  fim  and  organic/mineral 
grains  of  radii  0.03  ^m  have  Poynting-Robertson  timescales  of  3000  yr  on 
the  average,  whereas  graphite  particles  have  considerably  shorter  lifetimes. 
The  longest  Poynting-Robertson  timescale  is  calculated  for  porous  dielectric 
grains,  with  t  >  104  yr.  A  zodiacal  cloud  resulting  from  a  cometary  evapora¬ 
tion  episode  would  thus  have  an  average  lifetime  of  ~  104  yr. 

Modelling  this  temporary  zodiacal  cloud  as  a  disc  of  mass  5  x  1021  g, 
radius  1  AU  and  thickness  0.2  AU,  it  is  found  that  ~  109  tons  of  dust 
are  swept  up  by  the  Earth  annually  over  a  few  millennia  when  the  comet 
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is  active.  A  porous  mineral/organic  IDP  particle  of  radius  0.1  -  0.3  gm  and 
density  1  g  cm~3  has  a  settling  time  through  the  atmosphere  of  about  3-10  yr 
and  such  grains  are  efficient  scatterers  in  the  optical  wavelength  range.  It  is 
readily  estimated  that  the  Earth  gathers  up  a  stratospheric  dust  veil  whose 
optical  depth  may  fluctuate  in  the  range  0.05  ;$  r  <  3  during  the  cometary 
active  lifetime  of  103-104  yr. 

For  a  mixture  of  micron  and  submicron  grains  incident  on  the  Earth’s 
upper  atmosphere,  the  submicron  component  would  tend  to  be  optically 
dominant.  The  mean  scattering  angle  of  a  0.1  gm  particle  is  ~  90°  [?],  almost 
independently  of  composition,  and  its  settling  time  through  the  stratosphere 
is  about  a  decade  [26].  Because  the  settling  time  of  /im  grains  is  about  1 
year,  it  follows  that  if  equal  masses  of  gm  and  0.1  gm  grains  are  incident  on 
the  upper  atmosphere,  the  smaller  sizes  will  accumulate  in  the  stratosphere 
relative  to  the  larger  ones.  Further,  very  fragile  sub-cometary  meteoroids  of 
~metre  dimensions  may  disintegrate  to  ~  0.1  gm  particles  during  their  pre- 
atmospheric  descent  [27],  and  may  compete  with  the  preformed  zodiacal  dust 
as  a  major  source  of  optical  depth.  One  thus  expects  that,  in  the  presence  of 
a  disintegrating  giant  comet,  the  Earth  will  be  enveloped  in  a  highly  reflective 
dust  cloud. 


Glaciation  and  Greenhouse 

Without  the  greenhouse  effect  the  Earth’s  mean  temperature,  averaged 
over  latitude,  day  and  night,  and  land  and  sea,  is  given  by  [F©(1  —  A)/ac]1/4, 
where  Fq  =  1.37  x  10-6  erg  cm~2  s*1  represents  the  solar  energy  flux  at  the 
Earth’s  distance  from  the  Sun.  Here  A  is  an  averaged  value  for  the  Earth’s 
albedo,  c  is  the  speed  of  light  and  a  =  7.565  x  10 ~15  erg  cm~3  deg~4  is  the 
radiation  density  constant.  Thus  for  an  albedo  of  0.4  the  Earth  would  have 
a  mean  temperature  ~245  K  or  -28  C.  Opacity  sources  are  highly  wavelength 
dependent,  and  we  seek  an  approximation  with  the  virtue  of  physical  recti¬ 
tude  rather  than  attempt  to  set  up  a  supposedly  accurate  computation  in 
which  approximations  of  uncertain  physical  validity  are  nevertheless  made  in 
the  end!  Thus  if  we  divide  the  re-radiated  energy  at  ground  and  sea  level 
into  two  essentially  equal  halves  about  a  central  wavelength  at  A  =  13.5 gm 
(say),  then  we  may  suppose  that  half  the  emission,  longward  of  13.5  gm,  is 
completely  blocked  by  the  heavy  opacity  of  greenhouse  gases  while  the  re- 
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mainder,  shortward  of  13.5  /J.m,  is  completely  free  to  escape.  It  follows  that 
the  greenhouse  effect  raises  the  Earth’s  mean  temperature  by  a  factor  of  21/4, 
i.e.  to  about  292  K  (19  C)  rather  than  245  K,  a  result  agreeing  closely  with 
experience. 

In  practice,  water  vapour  and  CO2  are  the  main  greenhouse  gases.  The 
CO2  produces  essentially  the  whole  of  its  effect  through  absorption  in  the 
infrared  over  13.5//m  to  17.5/xm.  Because  the  blocking  by  carbon  dioxide 
over  this  interval  is  large,  the  band  having  steeply-falling  wings,  additions 
of  CO2  have  only  a  second-order  influence  on  the  greenhouse  effect  and  are 
inconsequential  compared  to  the  major  factor,  water  vapour,  controlling  the 
Earth’s  climate.  There  is  weak  absorption  by  water  vapour  shortward  of 
A  =  13.5  fim  while  its  blocking  effect  rises  all  the  way  from  17.5  to  ~  lOO^m. 
In  other  words,  there  is  a  partial  window  at  17.5  ;$  A  ;$  20  fim,  which  roughly 
compensates  the  blocking  shortward  of  A  =  13.5 fim,  essentially  justifying 
the  adopted  approximation.  Evidently  then,  it  is  water  vapour,  that  of  a 
‘standard  atmosphere’  implying  1  cm  cm-2  of  precipitable  water,  which  raises 
the  Earth’s  temperature  by  ~40  C  and  without  it,  the  Earth  would  now  be 
locked  in  a  permanent  ice  age  [28]. 

Reducing  the  water  content  of  the  atmosphere  to  a  few  mm  cm-2  weakens 
the  greenhouse,  dropping  the  Earth’s  mean  temperature  (for  the  same  A)  to 
~280  K,  which  corresponds  closely  to  what  is  required  for  typical  ice  age 
conditions.  It  follows  that  a  reduction  of  the  average  water  content  by  about 
two  thirds,  while  maintaining  the  albedo,  would  produce  an  ice  age.  The  same 
result  would  be  achieved  if  A  increased  from  0.4  to  0.5  or  if  the  insolation 
of  the  upper  atmosphere  declined  by  about  16  percent.  Such  albedo  and 
insolation  are  readily  attained  when  the  Earth  intersects  material  in  orbit 
associated  with  the  disintegration  of  a  giant,  active  comet.  Very  large  comets, 
it  appears,  will  trigger  the  onset  of  glaciations. 

Ice-age  conditions  have  generally  been  dry  and  cold.  The  great  deposits 
of  loess  (wind-blown  soil)  in  Eastern  Europe  and  China  imply  a  climate  that 
was  dusty  in  the  lower  atmosphere,  implying  a  low  precipitation  rate.  Low 
precipitation  is  not  a  handicap  in  the  accumulation  of  large  glaciers,  which 
will  grow  even  at  annual  precipitation  rates  as  little  as  a  few  cm/yr,  provided 
that  the  temperature  is  low  enough  to  prevent  summer  melting. 

During  the  ice  ages  the  whole  Earth  was  cooled,  including  the  tropics. 
This  is  proved  by  glaciers  extending  down  to  about  10,000  feet  on  tropical 
mountains  which  do  not  at  present  hold  glaciers,  such  as  those  on  Hawaii. 
The  need  for  the  whole  Earth  to  be  appreciably  cooled  is  difficult  to  rec- 
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oncile  with  ice-age  theories  depending  solely  on  small  oscillations  of  the 
Earth’s  rotation  axis  relative  to  the  ecliptic  plane,  and  small  oscillations  in 
the  Earth’s  orbital  eccentricity.  Neither  of  these  effects  produces  any  change 
in  the  amount  of  solar  energy  incident  on  the  Earth  and  so  could  not  lead 
to  widespread  cooling.  Oscillations  of  tilt  merely  produce  slight  latitudinal 
variations  in  the  incidence  of  solar  energy,  which  are  in  any  case  much  smaller 
than  the  transport  in  latitude  of  heat  by  atmospheric  storms  and  ocean  cur¬ 
rents.  The  transport  of  oceanic  heat  towards  the  poles  gives  a  far  larger 
effect  and  would  easily  buffer  slight  latitude  variations  of  insolation.  Oscilla¬ 
tions  in  eccentricity  of  the  Earth’s  orbit  produce  small  shifts  of  solar  energy 
between  one  geographical  hemisphere  and  the  other,  and  so  should  tend  to 
cool  one  hemisphere  and  warm  the  other.  But  ice-ages  occur  contempora¬ 
neously  in  both  hemispheres,  not  alternately,  a  disproof  that  was  already 
understood  more  than  half  a  century  ago.  The  presence  of  the  Milankovitch 
cycles  is  currently  controversial:  the  most  accurate  climatic  chronology  cov¬ 
ering  the  last  500,000  yr,  obtained  from  calcite  deposits  in  the  Devil’s  Hole 
fissure  in  Nevada,  appears  to  indicate  a  chaotic  behaviour  for  the  climate 
over  this  period  [29].  However  the  marine  pelagic  record  appears  to  support 
the  100,000-year  cycle  [15];  this  may  be  due  to  a  combination  of  planetary 
resonances  resulting  in  periodic  giant  comet  dusting. 

If  we  were  to  imagine  a  cold,  dry  atmospheric  state  being  brought  about 
today,  evaporation  from  the  relatively  warm  surface  layers  of  the  ocean  would 
quickly  resupply  water  vapour  to  an  amount  of  1  cm  of  precipitable  water 
per  cm2,  and  the  cooling  due  to  a  reduced  greenhouse  effect  would  quickly  be 
gone.  Thus  it  is  the  heat  of  the  ocean  which  saves  us  from  the  possibility  of 
an  immediate  onset  of  ice-age  conditions.  Reckoning  the  heat  of  the  ocean 
as  being  the  energy  content  above  freezing  point,  which  can  be  thought  of  as 
available  heat,  almost  all  is  contained  in  a  surface  layer  with  depth  no  greater 
than  a  few  hundred  metres,  the  amount  being  equivalent  to  a  supply  of 
sunlight  over  a  time  interval  of  a  few  years,  say  3  to  5  years.  It  is  because  the 
ocean  has  this  back  storage  of  heat  that  we  do  not  drop  almost  immediately 
into  an  ice  age. 

In  distant  geological  periods  the  heat  storage  in  the  ocean  was  consid¬ 
erably  greater  than  it  is  at  present.  Today  the  ocean  waters  are  close  to 
freezing,  whereas  only  50  million  years  ago  the  bottom  temperature  was 
~  15°C  and  the  available  oceanic  heat  was  then  equivalent  to  a  50  year  sup¬ 
ply  of  sunlight.  The  difference  has  been  attributed  to  drifting  continents, 
especially  by  the  positioning  of  Antarctica  and  Greenland  at  or  close  to  the 
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poles.  Melt  water  from  Arctic  glaciers  has  gradually  filled  the  lower  ocean 
with  water  close  to  freezing,  greatly  reducing  the  margin  of  safety  against 
ice-age  conditions  developing.  The  past  two  million  years  have  essentially 
been  a  continuing  ice-age  therefore,  broken  occasionally  by  short-lived  in¬ 
terglacials.  But  clearly,  from  this  perspective,  the  prospect  of  a  drift  away 
from  an  ice  age,  through  an  enhanced  greenhouse  effect,  may  be  of  less  con¬ 
cern  than  that  of  a  drift  back  into  an  ice  age.  Thus  we  need  sustain  the 
insolation  at  a  significantly  lower  value  than  at  present  for  only  several  years 
to  lock  enough  water  vapour  into  ice  to  create  the  permanently  cold,  dry 
atmospheric  conditions  of  an  ice  age. 


Changes  in  the  Earth’s  albedo 

A  remarkable  feature  of  the  Earth’s  albedo  is  that  it  may  very  easily 
be  raised  to  values  close  to  unity.  The  mass  extinction  coefficient,  through 
the  scattering  back  into  space  of  sunlight,  produced  by  dielectric  crystals 
with  radii  of  a  few  tenths  of  a  micrometer,  is  ~  3000  cm2  g~x  (backscattering 
amounts  to  about  10%  of  the  total  scattering  [40]).Such  an  albedo  change 
would  arise,  for  example,  if  even  a  very  small  fraction  of  even  a  very  dry 
atmosphere  were  to  condense  into  ice  crystals.  Thus  a  crystallization  of  only 
0.1  percent  of  the  water  in  a  very  dry  atmosphere  (say  with  only  1  mm  of 
precipitable  water)  would  yet  contribute  about  0.3  to  A.  Essentially  no  water 
must  be  condensed  into  ice  crystals  if  A  is  to  be  appreciably  less  than  unity; 
otherwise  the  Earth  would  appear  from  the  outside  as  an  intensely  bright 
white  planet  with  an  albedo  even  higher  than  that  of  Venus,  while  below  the 
haze  of  ice  crystals  it  would  be  exceedingly  cold  at  ground  level.  The  same 
result  would  arise  from  the  injection  of  ~  1014  g  of  ~  0.1  gm  dielectric  parti¬ 
cles  to  the  upper  atmosphere.  This  might  easily  happen  if  the  zodiacal  cloud 
were  to  be  flooded  with  the  debris  of  an  exceptionally  large,  disintegrating 
comet  in  a  short-period  orbit,  or  if  asteroidal  debris  from  such  a  comet  were 
to  strike  the  Earth.  We  discuss  here  the  ice  and  dust  mechanisms  in  turn. 

The  saving  grace  in  the  case  of  ice  crystals  is  that  they  do  not  form  in 
supersaturated  water  vapour  except  at  very  low  temperatures,  say  -50  C  [30]. 
For  the  Earth’s  emission  into  space  of  radiation  at  wavelengths  longer  than 
A =20  gm  we  can  think  of  a  photosphere  at  which  the  optical  depth  out  into 
space  is  of  order  unity.  If  only  radiation  were  involved  in  determining  the 
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water  vapour  temperature  at  this  photosphere  the  temperature  would  be  of 
order  290 r-1/4,  where  r  was  the  optical  depth  from  ground  level  up  to  the 
photosphere,  suitably  averaged  at  wavelengths  >  20  \im.  In  a  typical  atmo¬ 
sphere  r  would  be  about  10,  leading  to  a  photospheric  temperature  for  water 
vapour  (and  hence  for  surrounding  air)  of  as  little  as  163  K,  i.e.  -110  C,  far 
below  that  needed  for  ice  crystal  formation.  The  circumstance  that  ice  crys¬ 
tals  do  not  form  profusely  except  under  special  circumstances  in  Antarctica 
shows  that  calculating  for  radiation  alone  cannot  be  correct.  A  convective 
transport  of  energy  from  ground-level  to  the  water-vapour  photosphere  is 
required.  This  cannot  be  carried  by  air  movements  but  must  come  from  the 
upward  transport  of  the  water  vapour  itself.  To  keep  the  photospheric  water 
vapour  temperature  above  -50  C,  and  so  to  prevent  ice  crystal  formation,  the 
transport  of  water  vapour  must  be  such  as  would  lead  to  an  annual  precipi¬ 
tation  rate  of  about  50  cm  [30].  For  comparison,  the  present-day  worldwide 
average  of  the  precipitation  rate  is  about  80  cm  of  rain,  sufficient  to  prevent 
ice  crystal  formation,  but  not  by  a  wide  margin.  Let  the  world  climate  de¬ 
cline,  however,  sufficiently  for  the  surface  layers  of  the  ocean  to  cool  to  the 
point  where  an  annual  average  rainfall  of  50  cm  cannot  be  maintained,  and 
the  consequent  formation  of  an  atmospheric  haze  of  ice  crystals  would  plunge 
the  Earth  immediately  back  into  an  ice-age.  Such  sensitivity  of  the  Earth’s 
albedo  due  to  ice  crystallization  would  appear  to  indicate  an  additional  ten¬ 
dency  on  the  part  of  the  atmosphere-ocean  circulation  to  produce  rapid  ice 
age  conditions  in  the  presence  of  externally  injected  dust  [31]. 

The  KT  Event  Revisited 

Amongst  the  effects  which  were  predicted  [32]  as  a  consequence  of  gi¬ 
ant  comet  incursions  into  short-period,  Earth-crossing  orbits  were  multiple 
impacts,  high  concentrations  of  extraterrestrial  material,  prolonged  climatic 
deterioration,  ocean  regressions  and  a  complex  depositional  history.  These 
effects  are  over  and  above  the  prompt  effects  associated  with  large  impacts, 
and  permit  a  discrimination  between  the  stray  impact  and  cometary  hy¬ 
potheses.  Since  the  KT  boundary  has  been  intensively  studied  for  the  past 
decade,  it  provides  a  good  test  between  them.  Thus  amino  acids  of  probable 
extraterrestrial  origin  appear  to  have  been  laid  down  over  a  ~0.1  Myr  pe¬ 
riod  around  the  boundary.  Any  such  molecules  would  have  been  destroyed 
in  the  ~  105K  fireball  [33,  34].  The  data,  however,  are  consistent  with  a 
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protracted  input  of  extraterrestrial  organics  in  the  form  of  submicron-sized 
IDPs.  Another  observation  suggesting  a  prolonged  period  of  dust  input  is 
that  the  probable  extraterrestrial  component  of  minerals  at  several  impact 
sites  appears  to  be  one  or  two  powers  of  ten  too  high  to  be  consistent  with 
the  dilution  factor  ~  103  expected  for  impacts  (e.g.  ~100%  for  a  basal  layer 
of  clay  at  Woodside  Creek,  New  Zealand:  ref.  [35]).  Likewise  the  mass  of 
iridium  deposited  worldwide  is  overabundant  by  a  similar  factor  in  relation 
to  the  probable  size  of  the  Chicxulub  impactor  [36]. 

As  we  have  already  seen,  porous  submicron  grains  ( a  ~  0.1/xm)  of  either 
mineral  or  organic  composition  expelled  from  comets  are  mostly  retained  in 
bound  orbits  and  have  Poynting- Robertson  lifetimes  of  ~0.1  Myr  at  an  or¬ 
bital  radius  of  ~1  AU.  Such  particles  injected  into  the  zodiacal  cloud  will 
be  accreted  by  the  Earth  over  the  timescale  as  witnessed  in  the  distribution 
of  extraterrestrial  amino  acids  (AIB).  On  the  other  hand,  iridium-bearing 
graphite/metallic  grains  of  sizes  0.01  fim  (appropriate  for  supernova  con¬ 
densates)  will  be  expelled  from  the  inner  solar  system  on  a  much  shorter 
timescale.  On  this  basis  it  is  possible  to  understand  why  the  iridium  peak 
at  the  K/T  boundary  clays  is  considerably  sharper  than  the  0.1  Myr- wide 
profile  of  the  extraterrestrial  amino  acids  [37]. 

A  shroud  of  reflective  dust  accumulated  by  the  Earth  and  maintained  for 
some  0.1  Myr  would  lead  to  oscillations  of  climate  and  environmental  stresses 
causing  an  extended  episode  of  extinctions  of  species.  Such  a  picture  is  con¬ 
sistent  with  terrestrial  mean  annual  temperatures  estimated  from  studies  of 
fossilised  leaves.  There  is  an  indication  of  a  fourfold  rise  in  precipitation  and 
10°  increase  in  temperature  for  a  period  of  0.5-1  Myr  around  the  K/T  bound¬ 
ary  [41].  A  wide  range  of  geological  and  paleontological  data  seems  to  require 
the  combination  of  a  fairly  sharp  extinction  spurt  centred  at  65  Myr  as  well 
as  a  more  extended  episode  of  stepwise  extinctions.  For  dinosaurs,  however, 
a  consensus  is  emerging  that  there  was  no  gradual  decline  in  diversity  of 
genera  towards  the  end  of  the  Cretaceous,  but  rather  a  sudden  extinction 
consistent  with  a  concentrated  cluster  of  bolide  impacts  [42].  On  the  other 
hand,  groups  such  as  the  rudists  and  inoceradists  (bivalves)  seem  to  have 
disappeared  ~1  Myr  before  the  mass  extinctions  on  land.  Similarly,  Creta¬ 
ceous  foraminifera  and  dynocysts  in  K/T  boundary  clays  start  disappearing 
well  before  the  Ir  enhancement  begins  [43].  Depositional  and  palynological 
evidence  of  a  double  impact  layer  in  the  western  USA  has  been  presented 
by  a  number  of  workers  [38].  The  intervals  given  vary  from  a  few  months 
to  ~100  yr,  which  are  several  orders  less  than  those  expected  for  a  comet 
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shower,  but  readily  expected  from  a  dense  swarm  of  cometary  debris  within 
the  debris  stream  of  a  short-period,  Earth-crossing  orbit.  Multiple  bolide 
strikes  are  also  indicated  by  the  extensive  distribution  of  soot  and  products 
of  resinous  combustion  found  in  the  K/T  boundary  clay  layer  [44,  45].  The 
evidence  that  a  quarter  of  the  entire  biomass  was  combusted  is  inconsistent 
with  a  single  impact,  but  favours  the  idea  of  multiple  bolide  impacts  leading 
to  extensive  forest  fires. 


Effects  on  Civilization  and  Culture 

During  the  last  twenty  years  or  so,  there  has  grown  up  the  idea  that 
random  impacts  due  to  stray  asteroids  are  the  dominant  external  influence 
on  evolution  [11].  However  if,  rather,  such  dramatic  events  as  mass  extinc¬ 
tions  of  species  and  global  climatic  catastrophes  are  attributed  to  swathes 
of  cometary  debris  a  la  Shoemaker-Levy,  it  would  be  naive  to  gloss  over  the 
possible  implications  of  the  same  process  in  relation  to  human  culture.  The 
Tunguska  event  of  1908  seems  most  likely  to  have  been  caused  by  a  bolide 
of  ~100  m  diameter  exploding  at  a  height  of  some  8  km.  The  resulting 
blast  wave  felled  trees  over  a  distance  of  40  km,  charring  them  for  up  to 
15  km  from  the  centre  of  impact.  Estimates  of  a  300-100  year  timescale  for 
successive  Tunguska- type  collisions  are  based  on  lunar  cratering  data  [46]; 
this  calculated  rate  has  already  triggered  interest  in  projects  such  as  Space- 
guard.  However  the  lunar  maria  are  ~3.9  Gyr  old,  and  there  is  no  obvious 
reason  why  the  contemporary  impact  rate  should  bear  much  relation  to  that 
averaged  over  this  long  time  interval.  Recent  satellite  observations  of  the 
Earth  between  1975  and  1992  have  revealed  that  some  136  sub-Tunguska 
bolides  impacted  the  Earth’s  upper  atmosphere  during  this  interval,  yielding 
a  rate  of  a  few  per  annum  for  objects  10-30  m  in  diameter.  For  100  m  sized 
Tunguska-like  bolides  the  current  impact  rate  could  well  be  one  in  30-100 
years  [47,  34]. 

One  Tunguska-like  strike  per  century,  despite  its  attendant  horrors,  would 
have  little  sociological  impact.  But  if  in  times  past  similar  strikes  occurred 
at  the  rate  of  several  tens  per  annum  the  effect  upon  our  social  systems 
would  unquestionably  have  been  profound.  Such  collisions  are  possible  if 
the  Earth  intercepts  a  debris  stream  from  a  disintegrated  giant  comet,  as  we 
have  seen.  A  wealth  of  historical  data  exists  [34,  3]  to  support  the  hypothesis 
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that  fragmentation  of  a  particular  comet  began  some  20,000  years  ago  and 
the  interception  of  its  fragments  on  a  periodic  basis  has  led  to  events  that 
moulded  our  religions,  mythology,  beliefs  and  history.  The  beginning  of  the 
present  interglacial  period  is  marked  by  a  very  sharp  rise  in  temperature  at 
~13,000  BP,  followed  by  cooling,  and  a  further  sharp  rise  at  ~10,500  BP 
which  is  subsequently  maintained.  If  this  cooling  event  was  occasioned  by  an 
intense  spurt  of  cometary  bolide  impacts  due  to  fragments  of  a  giant  comet  in 
Earth-crossing  orbits,  similar  collision  episodes  with  declining  intensity  may 
have  continued  repeatedly  throughout  history.  One  could  regard  otherwise 
enigmatic  events  in  history  such  as  the  sudden  collapse  of  the  Indus  Valley 
Civilization  of  Mohendojaro  and  of  the  Old  Kingdom  in  Egypt  (accompanied 
by  the  most  puzzling  phenomenon  of  Pyramid  building),  both  occurring  at 
~2500  BC,  as  fitting  well  with  precession  of  the  primary  orbital  nodes,  the 
countdown  to  intersection  ca  500  ±  2500n  AD,  n=0,l,2...  [3]  and  the  (Tau- 
rid)  cometary  collision  picture  in  general.  Such  intersections  accompanied  by 
widespread  global  cooling  are  predictably  complemented  by  extended  periods 
of  global  warming  ca  1750  ±  2500n  AD  apparently  characterised  now  by  the 
recorded  ice-rafting  of  oceanic  sedimentary  debris  at  these  epochs  through¬ 
out  the  Holocene  [49].  The  oldest  celestial  myths,  which  involve  battles  for 
supremacy  between  gods  in  the  sky,  may  date  from  the  third  millennium  BC 
although  they  survived  at  least  through  to  Homer  and  Hesiod  ~800  BC  [34]. 
The  next  episode  of  violent  collisions  at  ~1000  BC  may  well  have  generated 
Old  Testament  accounts  such  as  the  destruction  of  Jericho.  On  this  basis  the 
most  recent  episode  of  severe  Tunguska-type  collisions  may  have  occurred  at 
~500  AD,  a  time  which  tallies  with  bizarre  phenomena  that  seem  to  have 
accompanied  the  end  of  the  Roman  Empire.  Gibbon  [48]  refers  to  a  ‘fever  of 
the  Earth  that  raged  with  uncommon  violence  during  the  reign  of  Justinian 
(AD  527-565)... Each  year  is  marked  by  the  repetition  of  Earthquakes,  of  such 
duration  and  severity  that  Constantinople  has  been  shaken  for  above  forty 
days...’.  As  the  Roman  Empire  collapsed,  so  did  that  of  the  Guptas  in  India, 
which  it  seems  was  torn  asunder  by  the  revolt  of  the  Huns.  Further,  W.M. 
Smart  refers  to  Islamic  text  at  a  similar  time  which  states:  ‘In  the  yeax  599 
on  the  last  day  of  Moharrem,  stars  shot  hither  and  thither  and  flew  against 
each  other  like  a  swarm  of  locusts;  people  were  thrown  into  consternation 
and  made  supplication  to  the  Most  High.’ 

Purely  scientific  evidence  from  dendrochronology  has  also  been  adduced 
to  support  the  idea  of  an  externally  caused  ecodisaster  occurring  at  around 
540  AD  [50].  From  studies  of  tree  ring  thicknesses  corresponding  to  the 
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early  decades  of  the  6th  century  AD,  it  has  been  found  that  a  major  dip  in 
the  Earth’s  temperature  occurred  over  the  entire  period  AD  536-546.  The 
competing  idea  that  a  volcanic  eruption  was  responsible  for  a  dust  shroud 
that  lowered  the  mean  temperature  and  reduced  tree  growth  for  a  decade  or 
more  does  not  accord  with  the  lack  of  an  acid  signal  in  Greenland  ice-drills 
encompassing  the  same  age.  Furthermore,  volcanic  dust  is  known  to  settle  in 
a  couple  of  years  at  the  outset,  so  cannot  easily  explain  the  protracted  episode 
of  cooling  that  has  been  found.  We  note  that  contemporary  literature  also 
concurs  with  dendrochronological  data.  It  has  been  stated  that  ‘the  sun  was 
dark  and  its  darkness  lasted  eighteen  months. ..the  sun  appears  to  have  lost 
its  wonted  light  and  appears  of  a  bluish  colour.. .fruits  did  not  ripen. ..cold  and 
drought  finally  succeeded  in  killing  off  the  crops  in  Italy  and  Mesopotamia 
and  led  to  terrible  famine  in  the  following  years.’  The  scene  was  surely  set 
for  widespread  mayhem  and  the  collapse  of  empires  that  followed. 
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Table  1.  Periodicity  in  the  terrestrial  record. 

Left-hand  column:  peaks  from  a  power  spectrum  analysis  of  all  global  events  (excluding 
those  more  recent  than  8  Myr  BP).  Center  column:  individual  peaks  in  the  mass  extinction 
record.  Right-hand  column :  peaks  from  a  PSA  of  all  global  geological  events  (excluding 
<  8  Myr  BP). 


‘All  data’  peaks 

Extinction  peaks 

Geological  peaks 
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11.9 

11.2 

10.3 

37.8 

36.6 

36.6 

63.7 
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62.9 

88.8 
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89.3 
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141.4 

144.0 

141.9 

167.3 

176.0 

168.2 

193.2 

193.0 

194.5 

219.0 

216.0 

220.9 

245.0 

245.0 

247.2 
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Figure  captions 

Fig.  1.  Top:  Percentage  of  marine  genera  extinctions  from  mid-Permian 
to  Recent,  according  to  Raup  &  Sepkoski  [1].  Excluding  the  recent  surge 
of  activity  attributed  to  Gould’s  Belt  disturbance,  the  best-fit  periodicity  is 
P  =  12.2  +  25.9n,  n  =  0, 1,2, ..  (P  in  millions  of  years  BP). 

Bottom:  Major  geological  events  [12],  plotted  with  exponential  smoothing 
within  a  4  million  year  window.  Best-fit  periodicity,  excluding  recent  activity, 
of  P  =  10.3  +  26. 3n,  n  =  0, 1, 2, ...  is  marked  by  vertical  dashes. 

Fig.  2.  Stability  of  the  combined  extinction  and  geological  data  against 
progressive  truncation  of  the  data,  Tmax  the  upper  limit  (the  recent  out-of- 
phase  surge  is  again  excluded).  Left:  Variation  of  period  P  and  phase  <j>. 
Right:  Behaviour  of  typical  synthetic  data  simulating  the  real  data  (with 
regard  to  trend  and  edge  effects)  but  with  each  datum  perturbed  randomly 
in  the  range  ±24  Myr.  Trials  with  such  data  favour  periodicity  over  chance 
at  a  confidence  level  ~  1  —  10~3. 

Fig.  3.  Top:  Modulation  of  Oort  cloud  flux  into  the  planetary  system  caused 
by  numerous  mini-showers,  generated  by  encounters  with  passing  dark  bod¬ 
ies.  In  this  case  the  assumed  perturbers  are  assumed  to  be  Population  I  black 
dwarves  which  penetrate  the  dense  inner  core  of  the  Oort  cloud,  while  the 
solar  orbit  is  taken  to  have  vertical  amplitude  80  parsec  in  the  Galactic  disc. 
The  adiabatic  tide  due  to  the  cumulative  effect  of  remote  perturbers  yields  a 
similar  flux  amplitude.  Bottom:  As  above  but  with  4  million  year  smoothing 
corresponding  roughly  to  the  diffusion  time  of  comets  within  the  planetary 
system  [13]. 

Fig.  4.  Representation  of  Gould’s  Belt  and  related  gas,  adapted  from 
Olano,  C.A.,  Astron.  Astrophys.  112,  95  (1982).  The  system  is  expanding 
and  comprises  an  outer  ring  of  young  stars  which  has  detached  from  an  inner 
gas  ring.  Its  mass  M  106  A/©,  and  the  current  distance  to  its  center,  166  pc 
(Olano)  or  80±40  pc  (Comeron,  F.  &  Torra,  J.,  Astron.  Astrophys.  241, 
57  (1990)),  are  critical  in  respect  of  the  inferred  relative  amplitudes  of  the 
impulsive  and  steady  components  of  the  galacto-terrestrial  cycle.  The  solar 
system  appears  to  have  passed  through  the  rim  of  the  Belt  <  10  Myr  BP. 
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Fig.  5.  The  ratio  of  radiation  pressure  to  gravity  for  spherical  particles  near 
the  Sun.  The  curve  G  refers  to  graphite;  I  refers  to  iron;  D100,  D60,  D40 
refer  to  a  sequence  of  non-porous  grains,  60%  volume-filled  porous  grains  and 
40%  volume-filled  porous  grains  respectively.  The  dotted  and  dashed  lines 
show  critical  P/G  values  for  retention  in  bound  orbits. 

Fig.  6.  The  characteristic  Poynting-Robertson  lifetimes  for  particles  in 
heliocentric  orbit  at  1  AU.  The  symbols  marking  the  curves  refer  to  the  same 
cases  as  in  Figure  5. 
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Appendix  A: 


the  nature  of  the  cometary  hazard 


The  celestial  hazard  ostensibly  feared  most  by  our  ancestors  during  enhancements  of  the 
global  fireball  flux  was  a  sustained  bombardment  by  larger  cometary  debris.  Included  here 
is  the  Prologue  of  "The  Cosmic  Winter"  by  Victor  Clube  and  Bill  Napier  (Basil 
Blackwell),  Oxford  1990),  conceived  before  the  Cold  War  came  to  an  end  and  providing  a 
sketch  of  the  feared  conditions  translated  to  a  contemporary  environment.  This  sketch 
anticipated  the  tidal  break-up  of  Comet  P/Shoemaker-Levy  9  and  its  subsequent 
bombardment  of  Jupiter  by  four  and  six  years  respectively. 
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